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ABSTRACT 

A combination of the X-ray imaging and SZ measurements of clusters permits 
an indirect determination of the radial temperature profiles of intracluster gas, 
which requires no assumption about the dynamical properties and the equation 
of state for clusters. A comparison of such a derived gas temperature with that 
given by the X-ray spectral analysis constitutes an effective probe of properties 
of intracluster gas. Using the available data of 31 clusters in literature, we 
have performed the first comparison of the central gas temperatures provided 
by the two methods. The good agreement between these two independent 
temperature estimates suggests that the distribution of intracluster gas is 
essentially consistent with isothermality characterized by a mean polytropic 
index of 7 = 0.9 ±0.1. 

Subject headings: galaxies: clusters: general — intergalactic medium — X-rays: 
galaxies 
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1. Introduction 

While the X-ray spectroscopic measurement is a powerful and unique way nowadays to 
obtain the temperature of hot, diffuse gas contained within clusters of galaxies, it yields an 
emission-weighted temperature rather than the true temperature distribution. The lack of 
detailed information on the temperature profiles in clusters is probably the major source 
of uncertainties in the present determinations of the total masses and baryon fractions of 
clusters, which further hinders the dynamical properties of clusters from the cosmological 
applications such as the estimates of the cosmological parameters and the test of various 
models of structure formation in the universe. An independent, complementary method of 
measuring the distribution of intracluster gas is thus desirable. Silk & White (1978) were 
the first to suggest the utilization of the nonparametric reconstruction of the radial profiles 
of density [n e (r)] and temperature [T(r)] of intracluster gas by combining the X-ray and SZ 
measurements. The basic idea is to inverse the observed X-ray and SZ temperature surface 
brightness profiles of clusters which have the functional forms of n^r)!" 1 / 2 ^) and n e (r)T(r), 
respectively. And a simple combination of the two functions gives straightforwardly the 
gas density n e (r) and temperature T(r). This method requires no assumption about the 
dynamical properties of clusters and the equation of state for intracluster gas, and can 
therefore be regarded as an ideal and ultimate tool to probe the gas distribution in clusters 
under spherical approximation. 

Despite its elegant mathematical treatment further developed by Yoshikaya & Suto 
(1999) based on some theoretical models, the pioneering suggestion of Silk & White 
(1978) has not yet been put into practice to date. This has been primarily limited by 
the instrumental sensitivity and resolution of detecting the temperature variations of the 
cosmic background radiation (CBR) behind clusters. So far, the marginal detections of 
the radial SZ temperature distributions have been reported only for a few clusters (see 
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Birkinshaw 1999), which can hardly be used for the purpose of reconstructing the gas 
temperature profiles because of the sparse data points and the large associated uncertainties. 
Nevertheless, we notice that the central measured or estimated SZ effects have been available 
for a great number of clusters with the past two decades' efforts, and these central CBR 
temperature decrements in the Rayleigh- Jeans limit are usually used for setting constraints 
on the cosmological parameters (H , Qm and i n conjunction with the X-ray imaging 
and spectral measurements. Yet, another application of the central SZ data is instead to 
estimate the central temperature of intracluster gas when combined with the X-ray imaging 
observation, although this does not completely achieve the original goal of Silk & White 
(1978). It provides an indirect measurement of the X-ray temperature at cluster centers, 
which can be directly compared with the result given by the X-ray spectral analysis. 
Additionally, because both X-ray emission and SZ effect result from the gas distribution 
projected along the light of sight, a comparison of the central temperatures indicated by the 
two methods may also reveal valuable information about the radial temperature variations. 
Using the published X-ray and SZ data in literature, we present for the first time such a 
comparison in this paper and discuss its implications for dynamical properties of clusters. 
Throughout the paper we assume H = 50 km s _1 Mpc^ 1 and Q = 1- 



The X-ray imaging observation provides a reliable measurement of the X-ray 
azimuthally-averaged surface brightness profile of a cluster which is usually described by 
the conventional f3 model (Cavaliere & Fusco-Femiano 1976) 

/ ^2 \ -3/3+1/2 



2. 



The model 




(1) 
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In the scenario of an optically-thin, thermal bremsstrahlung emission, the above form of 
S x (r) indicates (Cowie, Henriksen & Mushotzky 1987) 

, / r 2\ -3/3/2 

n e {r)T'l\r) = n e0 T 1/4 1 + -\ , (2) 



where n e and T are the electron number density and temperature, respectively. If we 
assume an equation of state, T(r) = T [n e (r)/n e o] 7_1 , we can write the electron number 
density as n e (r) = n e0 (l + r 2 /r 2 )~ s , in which 5 = 6/5/(3 + 7). The central electron number 
density n e0 , temperature T and the X-ray surface brightness So are connected by 

2 = ( 3m e hc 2 \ ( 3m e c 2 \ 1/2 An 1 ' 2 T(3(3) S (l + zf 
He0 { W )\27ckT ) ^r(3/3-l/2) r c ' U 

where /i e — 2/(1 + X), X — 0.768 is the primordial hydrogen mass fraction, g as 1.2 is the 
average Gaunt factor, T is the gamma function, and z is the cluster redshift. 

On the other hand, the CBR temperature decrement /increment at the cluster center 
predicted by the SZ effect is 

T CBR ^0 m e c 2 

x 2 e x ( , , x 



£(x) = — 4-rrcoth- , (5) 

^ V J 2(e x - l) 2 V 2) y ' 

where T CB r is the temperature of the present CBR, and x = hu/kT CBR is the dimensionless 
frequency. A straightforward computation using eqs.(2)-(5) as well as the polytropic 
equation of state gives 



kT \ 3/2 



(AT SZ \ 2 ( IW/2) \ 2 T(3/3-l/2) 
VTcbeJ \T{W'/2-l/2)) T(3/3) 

(6) 



m e c t 



87r 2 v / 3(3 / 3-3/2)5 2 <T T (l+2) 4 5 r c 

in which /3' = 3f3 r y/{3 + 7), and a is the fine structure constant. The above expression is 
often used in the determination of the Hubble constant for an isothermal gas distribution 
if we write r c = cIa(z)9 c where cIa is the angular diameter distance to the cluster at z. On 
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the other hand, given the Hubble constant and the polytropic index 7, all the parameters 
at the right-hand of eq.(6) are measurable from the X-ray imaging and SZ measurements, 
which will instead allow us to derive the central temperature of the X-ray emitting gas. 

3. Application to X-ray clusters 

Our cluster sample consists primarily of 21 nearby clusters studied extensively by 
Mason & Myers (2000) plus 10 high-redshift clusters with recent, reliable SZ data in 
the literature. Inclusion of a cluster is also subject to whether or not the good spectral 
temperature data are available in order to facilitate our comparison between different 
methods. We take the best-fit values of (3, r c and So (or L x ) as well as the emission-weighted 
temperature data directly from the literature. For the latter, we adopt the temperature 
given by a single-phase (s.p.) model and the value by excluding the cooling flows (c.l.) 
(White 2000), respectively. Since the X-ray spectrum is always observed over a finite 
bandpass, we assume an optically-thin, thermal bremsstrahlung emission model with the 
metal abundances of 30% solar to convert the observed So (or L x ) over the corresponding 
energy band into the bolometric value. We use the central SZ decrement /increment of AT SZ 
after the correction of the finite beam in each observation. All the input data of the X-ray 
and SZ measurements are summarized in Table 1. 

For an isothermal gas distribution, the jointly determined central temperature, T est , 
from the X-ray imaging and SZ measurements (eq.[6]) can be directly compared with the 
X-ray spectral result, T spec . We display in Fig.l our derived T est versus the observed T spec 
(s.p. model) for the 31 clusters in Table 1. A glimpse of Fig.l reveals that the two sets of 
data are roughly consistent with each other. The best-fit T spec -T est relation reads 

T est = io-o- 22 ^^ Ts i p 22 ± o.20 ) s p model . (7) 
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Tset = 10-0.36±0.26 T W9±0. 2Bj c J corrected) (g) 

in which the error bars are the 90% confidence limits which have taken the measurement 
uncertainties into account. We have also tested the dependence of the best-fit T spec -T est 
relation on the Hubble constant by adopting a larger value of Ho = 75 km s -1 Mpc -1 . This 
only leads to a minor modification to the above results. For instance, the best-fit relation 
for the s.p. model now reads T est = io-o.34±o.ii r ^±o.i2_ 

EDITOR: PLACE FIGURE |T| HERE. 

If the intracluster gas deviates from isothermality, our comparison of T est and T spec 
becomes to be a little complex because we need to recover the central temperature from 
the emission- weighted spectral fit. On the other hand, we can place a useful constraint 
on the polytropic index 7 if we demand that our derived central temperature from eq.(6) 
identify the spectral value after a proper correction to the emission-weighted temperature 
T predicted by 

^ oc T~ 1/2 g ff (T, v) exp(-hv/kT) J n 2 e dV, (9) 

where gff is the Gaunt factor of the free-free emission, v is the X-ray observing frequency, 
and the integration is performed over the whole cluster volume if we neglect the finite 
detection aperture. This compares with the theoretically expected spectrum for the hot 
intracluster gas 

— oc j n 2 e T~ 1/2 g ff {T,v)exp(-hv/kT)dV. (10) 

For each cluster with the reported T = T spec , we use the Monte Carlo simulations over the 
typical energy band 0.5 - 10 keV to obtain the central temperature T and the corresponding 
7 with the restriction of eq.(6). The resultant 7 value for each cluster is listed in Table 1, 
and a combined analysis of the results among the 31 clusters yields 

(7) = 0.91 ±0.14, s.p. model; (11) 
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( 7 ) = 0.93 ± 0.14, c.l. corrected, 



(12) 



which indicates that clusters are essentially consistent with isothermality at 90% significance 
level. 



Our joint analysis of the indirect and direct temperature measurements from the X-ray 
and SZ observations provides a useful clue to resolving the temperature profile discrepancy 
raised in recent years particularly by Markevitch et al. (1998), Irwin, Bregman & Evrard 
(1999), White (2000) and Irwin & Bregman (2000). Markevitch et al. (1998) claimed a 
significant temperature decline with increasing radius represented by a polytropic index of 
1.2-1.3, while the other investigators have essentially detected an isothermal temperature 
profile. Our present study gives an independent yet convincing evidence for the assumption 
of isothermality. The advantage of our algorithm over the previous analysis is that we 
have obtained the central temperatures of clusters without utilizing the X-ray spectral 
data, while a combination of both has allowed us to set a valuable constraint on the radial 
temperature profiles of the hot X-ray emitting gas. 

A joint determination of the central temperature of intracluster gas from the X-ray 
imaging and SZ measurements can be used at present as an independent and complementary 
method to examine the reliability and accuracy of the temperature measurement provided 
by the conventional X-ray spectral analysis, although our ultimate goal is to obtain the 
detailed radial profiles of electron number density and temperature when the SZ surface 
brightness can be precisely mapped (Silk & White, 1978; Yoshikawa & Suto 1999). Yet, 
we should point out that it may not be easy to identify the temperature structure of the 
hot gas from such a combined analysis, in the sense that the dominant component in both 
S x and AT SZ is characterized by the sharp decrease of the electron number density with 
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Discussion and conclusions 
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outward radius, while it is unlikely that the temperature profiles can demonstrate a similar 
variation. An accurate measurement of S x and AT SZ at outermost radii will thus be needed 
if one wants to apply the present technique to the reconstruction of temperature profiles 
of intracluster gas. Alternatively, the present technique requires that the X-ray surface 
brightness is well represented by the (3 model. Presence of cooling flows in the central 
regions of many X-ray clusters may yield an underestimate of the (3 parameter and the core 
radius r c and an overestimate of the central surface brightness So, which will in turn affect 
our estimate of T est through eq.(6). Indeed, an examination of Fig.l reveals that all the 
massive cooling flow clusters with M > 600 M yr^ 1 (A478, A1835, A2204 and Zw3146) 
demonstrate rather large scatters around the T est = T spec relation. Therefore, the present 
method is most applicable to morphologically simple clusters without cooling flows. 

This work was supported by the National Science Foundation of China, under Grant 
19725311. 
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Fig. 1. — The derived central temperatures are plotted against the spectral results (single- 
phase model) for a sample of 31 clusters. Solid line shows the best- fit relation, and dotted 
line indicates T est = T spec . The high-redshift (z > 0.1) and low-redshift (z < 0.1) clusters 
are represented by the filled squares and triangles, respectively. The marked clusters are the 
massive cooling flow ones with M > 600 M yr _1 (White 2000). 
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a Single-phase model 
b Cooling- flow corrected 
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